Particle model for skyrmions in chiral magnets: dynamics, pinning and creep 
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Recently spin textures called skyrmions have been discovered in certain chiral magnetic materials without 
spatial inversion symmetry, and have attracted enormous attention due to their promising application in spin- 
tronics since only a low applied current is necessary to drive their motion. When a conduction electron moves 
around the skyrmion, its spin is fully polarized by the spin texture and acquires a quantized phase; thus, the 
skyrmion yields an emergent electrodynamics that in turn determines skyrmion motion and gives rise to a fi- 
nite Hall angle. While Skyrmions behave as particles, no particle level description of their dynamics exists 
to date. Such a model would have tremendous impact on understanding skyrmion dynamics by theoretical 
analysis and computational modeling. Here we derive the equation of motion from a microscopic continuum 
model and obtain the short-range interaction between skyrmions, and the interaction between skyrmions and 
defects. Skyrmions also experience a Magnus force perpendicular to their velocity due to the underlying emer- 
gent electromagnetic field. We validate the equation of motion by studying the depinning transition using both 
the particle and the continuum models. By using the particle description, we explain the recent experimental 
observations of the rotation of a skyrmion lattice in the presence of a temperature gradient. We also predict 
quantum and thermal creep motion of skyrmions in the pinning potential. 

PACS numbers: 75.1().Hk, 75.25.-j, 75.30.Kz, 72.25.-b 



I. INTRODUCTION 

Spin texture called skyrmion was predicted to exist in cer- 
tain magnetic materials. 1 1 -3 1 Skyrmion crystals have been 
recently observed in MnSi, Feo.5Coo.5Si, and other B20 tran- 
sition metal compounds with small angle neutron scattering, 
Lorentz force microscopy and spin-polarized scanning tunnel- 
ing microscopy. l4,-i9J These spin textures become more sta- 
ble in thin films IIS fflTI . and they crystallize into a triangular 
lattice similar to that found for vortices in type II supercon- 
ductors. The typical size of a skyrmion is about 10 nm and 
the corresponding lattice constant is about 100 nm. As more 
skyrmion crystals are discovered in new materials, it is ex- 
pected that this state of matter will turn out to be a general 
form of magnetic ordering, existing ubiquitously in magnets 
without inversion symmetry. 

A promising set of spintronics applications arises from the 
fact that skyrmions can be driven by a spin-polarized current 
as a result of the spin-transfer torque. The weak current re- 
quired to move a skyrmion from the pinning center is 4 to 5 
orders of magnitude smaller than the current required to move 
the well-studied magnetic domain walls. fT2Ul4| Therefore, 
skyrmions can be manipulated with much less energy dissi- 
pation. A theoretical framework for understanding skyrmion 
dynamics is then crucial for applications. Current descrip- 
tions are based on continuum models that are difficult to solve 
analytically, and can be computationally intensive. Because 
skyrmions appear to have particle-like properties, the deriva- 
tion of a paiticle-based equation of motion provides a func- 
tional form for interactions between skyrmions, skyrmion- 
defect interactions, and the role of terms such as the Mag- 
nus force. The derivation of effective equations of motion 
for other systems, such as vortices in type-II superconduc- 
tors lITSji . has been crucial for understanding pinning and vor- 
tex dynamics in the flux-flow regime. 

In this article we derive a concise particle-like equation of 



motion for skyrmions. The emergent electromagnetism in- 
duced by the Berry phase leads to an additional Magnus force 
that strongly suppresses the depinning current by deflecting 
skyrmions away from the pinning centers. By applying the 
deiived equation of motion to the study of skyrmion lattice ro- 
tation in the presence of a temperature gradient, we reproduce 
recent experimental results. We use the same equation of mo- 
tion to investigate the quantum and thermal creep motion of a 
skyrmion in a pinning potential. Finally, we validate the parti- 
cle model by computing depinning transitions and comparing 
against results obtained with the original continuum model. 



II. EQUATION OF MOTION 

We consider a thin film of a chiral magnet with 
Dzyaloshinskii-Moriya (DM) interaction which supports 
skyrmions. |[T]-(3] [16] [iTl The magnetic moments are de- 
scribed by a unit vector n(r). The coiTesponding action for 
the magnetic moments n can be wiitten as 



da r 
J J 



dtdt'd r 



n(r) - n(f') 



f- f 



dt9{, (1) 



h^dj d^rdtz^'i{^^d,- ^J-Vjz, 



(2) 



where z = \z} is the spin coherent state defined as n ■ cr\z} - \z}- 
cr is the vector of Pauli matrices and d is the film thickness. 
Here, y - litis) with a the lattice constant and s the total 
spin. The first term in 5 b describes the Berry phase for the 
precession of a spin at r = (x, y). In the presence of con- 
ducting electrons, the electrons become fully polarized by the 
local moments n in the large Hund's coupling limit, as de- 
picted in Fig. [T] The second term in 5 b is responsible for the 
Berry phase that the electron picks up when it moves around 
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FIG. 1. (color online) Schematic view of an electron (yellow sphere 
with arrow) passing through a skyrmion (colored arrows). The spin 
of the electron follows the spin texture of the skyrmion, giving rise to 
an emergent magnetic field that couples the electronic orbital motion. 
The black arrows (dots) are the spin projection in the x-y plane and 
the spins in the ferromagnetic state are along the z axis. 



the skyrmion. The term proportional to a accounts for the 
Gilbert damping. The spin Hamiltonian is 



(3) 



The first term is the exchange interaction, the second term 
is the DM interaction, which breaks spatial inversion sym- 
metry, and the last term is the Zeeman energy. The exter- 
nal magnetic field is perpendicular to the film. Systems gov- 
erned by Eq. ([3]l support a skyrmion phase in an intermedi- 
ate magnetic field 0.2D^/J^^ < H,, < O.SDV^ex- El The 
skyrmion is characterized by the topological charge density 
Q{r) = / dr^n ■ (5,-n x 5vn)/(47r) = +1. 

According to Eqs. ([T} and (|2]), the spin dynamics is gov- 
erned by the Landau-Lifshitz -Gilbert equation lfT8l420l 



d,n = — (J ■ V)n - yn X Heff + ad,n 
2e 



X n. 



(4) 



with the effective magnetic field Heff = d'H/dn. 

Equations ([TJ-Q describe the skyrmion dynamics as well 
as its deformations. Skyrmions can be treated as particles 
as long as deformations of their internal structure remain 
small. In other words, a particle-like description assumes that 
skyrmions have a rigid internal structure. Such rigidity is de- 
termined by the frequency of the normal modes associated 
with small fluctuations of the spin texture around the station- 
ary state. Here it is important to not that the Zeeman and the 
DM terms remove any continuous symmetry except for trans- 
lations. Therefore, the only Goldstone mode arises from trans- 
lations of the rigid skyrmion. Modes associated with the inter- 
nal skyrmion structure always have a finite frequency. This fi- 
nite frequency gap increases with magnetic field and provides 
a natural justification for treating skyrmions as particles. To 



treat skyrmions as point particles, we make the following two 
approximations. We assume a skyrmion density such that the 
overlap between different skyrmions is small. We also assume 
that the structure of moving skyrmions is the same as that in 
the static case. The internal structure of skyrmions becomes 
iiTelevant under these conditions, which are satisfied in the 
low velocity region for certain magnetic fields. |'2T1 We first 
derive the equation of motion for a skyrmion in the particle- 
level description [see Methods A for derivation]: 



47ta 



V; = Fm + Fi + y F,,(r,- - r,-) + V F,,(r,. - r,-), (5) 



j 



where v, is the skyrmion velocity, rj - rj^ - An[^ (ir^(5^n)^]"' 
and n - X, y. T] X \ depends weakly on J^ilD for J^ilD » a. 
Equation (|5]l is the main result of the present work. The term 
on the left-hand side accounts for the damping of skyrmion 
motion, which is produced by the underlying damping of the 
spin precession. F^ = Any^^z x v, is the Magnus force per 
unit length, which is perpendicular to the velocity. F^ = 
2;r/je"'z X J is the Lorentz force, that arises from the emergent 
quantized magnetic flux Oq = hc/e carried by the skyrmion 
in the presence of a finite cuiTent. Fj, is the pairwise interac- 
tion between two skyrmions and is the interaction between 
skyrmions and quenched disorder. It is clear from Eq. (|5]l, 
that the skyrmion does not have an intrinsic mass. For thin 
films, the skyrmions are straight in the direction perpendicu- 
lar to the film, and the forces in Eq. (|5]l are defined per unit 
length. A similar equation of motion was considered before in 
the context of vortices of type II superconductors. [15,1 How- 
ever, the Magnus force is negligibly small for superconducting 
vortices. [15| 

The action for the particle model of Eq. ^ can be written 

as 



S B.p - 



' B.p ■ 



An 



U(r) ■ 



Ana 



dtdt' 



r(f) - r(f') 



r 
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(6) 



(7) 



where U (r) is the potential per unit length produced by other 
skyrmions and pinning sites. A unique feature is that x and 
y are conjugate variables, which is a hallmark of the Berry 
phase, as given by 5^ in Eq. (|2|. 

One can treat the Magnus force F^ as originating from 
an effective transverse magnetic field, - Anc/{yq), that 
couples to a charged moving particle with charge q. This 
emergent magnetic field originates from the Berry phase. 
The Magnus force does not produce work, but it affects the 
skyrmion trajectory. As we will see later, skyrmions are eas- 
ily deflected by pinning centers because of the Magnus force. 
This effect explains the very weak pinning that has been ob- 
served in different experiments. 

We estimate the force using the typical parameters for 
MnSi, [22] a ~ 2.9 A, 7ex ~ 3 meV/a, D * 0.3 meV/a^, 
a X O.l and s » 1. At a velocity v = 1 m/s, we estimate the 
dissipative force per unit length to be Fdis = Anav/iyrj) 
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FIG. 2. (color online) Profile of Q(r) obtained from a numerical solu- 
tion of Eq. l |10[ ) for different values of the magnetic field. 



5 X 10" N/m; the Magnus force per unit length is Fm ~ 
5 X 10"^ N/m. Thus Fm » /^dis- The repulsive force per unit 
length between skyrmions for d ^ 20 nm is Fss ~ 10"^ N/m 
at a separation = 10 nm (see Fig. [3]below). The Lorentz 
force per unit length at a current density J - 10^ A/m^ is 
Fl^ Ay. 10"^ N/m. Since the depinning current for skyrmion 
lattice is of the order of 10^ A/m^, lfT2Ul4l we thus estimate 
the pinning fore per unit length as * 4 x 10"' N/m. 




FIG. 3. (color online) Force between two skyrmions as a function 
of the separation in two different magnetic fields. Symbols are 
obtained from a numerical solution of Eq. (|4]l and lines are fits to 
Kiir^l^), with ^ = D/ ^H„J^^. Inset: stationary configuration of 
two skyrmions at rj = 3.6/cx/O. The vectors denote the and n,, 
components and the component is represented by the color scale. 



A. Derivations 

Spins precess collectively when a rigid skyrmion moves 
with velocity v, n5(r - yt), and their evolution is governed 
by the equation of motion 



hy 

dfiis - — (J ■ V)nj - yiij X H/ + an^ x c?,nj, 
2e 



(8) 



where H, =11^+ H^. Hj is the magnetic field produced by 
other skyrmions and is the field produced by defects. The 
effective field Hq = S'HIda.s due to the skyrmion does not 
contribute to H, because x Hq = for a rigid skyrmion. If 
we first multiply both sides of Eq. (|8]l by xns (cross product) 
and then by -S^n, (dot product), we obtain 



a y 
~v = — 
?7 An 



Fm + Fz, + J dr^n^r' - r) ■ V,n, 



(r) 



(9) 



after integrating over the area around the skyrmion. Here 
is the field component perpendicular to iij and we have used 
that J dp-dxiis ■ dyUs - for a rigid skyrmion. The interaction 
potential between a skyrmion at r and an another skyrmion 
at r' is t/js(r' - r) = - J dr"^ns(r - r") ■ H.5(r' - r") and the 

corresponding force isFjj - J (ir"^V,nj(r-r")-Hj ^(r'-r"). 
The self-energy of the skyrmion in the presence of defects is 
E,(r - r') = - / dr"^n,(r - r") ■ Hrf(r' - r"), where Hd(r) = 
/ex(r)V^ns/2 - D(r)V x n, + B. The pinning force is then 
given by = J dr"^Vn,(r-r") ■ H<^,i(r' - r"). Thus, Eq. ^ 
reduces to Eq. Q if we replace the integral by the interaction 
force. 

To calculate the interaction between skyrmions and the in- 
teraction between skyrmions and defects, we need to know 
the structure of a single skyrmion. An isolated skyrmion is 
described by ns(r, (p) = sin 0^ + cos 0z in the polar coordinates 
(r, <p) with ^ and z being the unit vectors along the correspond- 
ing axises. 0(r) is determined by minimizing in Eq. (|3]l. 



rdt.0 - dr0 + cos(26') ■ 



sin(26') B 
— ^+^rsin(0)- 
2r 2 



1 =0, (10) 



with the boundary condition 6{r - 0) - n and 0(r — > -i-oo) - 0, 
where we have renormalized the distance r as r — » rjiJ^^ID), 
and p - 2HaJex/D^- The profile of 6(r) for different /3 is 
shown in Fig. |2] There are two length scales associated with 
a skyrmion. decreases linearly in r for r <sc 1, while the 
asymptotic solution far away from the center of the skyrmion, 
r — > oo, is ~ Ko(r/^) with a healing length ^ = v/2//3. Here 



K() and Ki below are the modified Bessel functions. The spin 
recovers exponentially to the fully polarized state due to the 
finite energy gap in the spectrum of the spin wave excitations 
that is induced by the external field. One may define the core 
region of the skyrmion as 0(r < Rc) < 7t/2. 

Since the interaction between two skyrmions is induced 
by the overlap between both spin textures, it must depend 
on the length ^. To calculate the interaction between two 
skyrmions numerically (see Appendix B), we initially pinned 
the two skyrmions at a fixed separation, rd, by freezing the 
spins within a radius r < J^x/D during the time evolution dic- 
tated by Eq. (|4]), and calculated the energy as a function rd- 
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minimize their self-energy. 5) The non-uniformity of electron 
density in real solids (size of defects) is of the order of the 
inter-atomic length, ~ 0.1 nm, that is much smaller than 
the typical skyrmion size. Thus, the interaction between de- 
fects and skyrmions is extremely weak. This is one of the 
reasons why the pinning of skyrmions is very weak. 



III. APPLICATIONS 

We apply the particle model of skyrmions to study the ro- 
tation of the skyrmion lattice in the presence of a temperature 
gradient, as well as the creep motion. We also compare the 
particle and continuum models by considering the depinning 
transition. 



FIG. 4. (color online) (a) Interaction force between a skyrmion and 
a defect for different strength and size of defect. Lines are results 
obtained from full numerical calculations and symbols are results 
obtained by assuming a rigid skyrmion structure. 



A. Rotation of skyrmion lattice 

Recent neutron scattering experiments have shown that the 
skyrmion lattice rotates in the presence of a temperature gra- 



The results shown in Fig.[3]indicate that the interaction decays 
exponentially and it is well described by Fss ~ Ki(r^/^). 

We next proceed to study the interaction between 
skyrmions and defects. The electronic density is not homoge- 
neous in real systems leading to an inhomogeneous exchange 
interaction 7ex produced by the double-exchange mechanism. 
We model the defects by the following profile of J^^^: 



Jexir) = Jq 



1 +^7dexp[-|r-rrf,,|/^rf] 



(11) 



where 7^ characterizes the strength of the defects and ,■ is 
the pinning center The characteristic size of the defects, ^d, 
is comparable to the inter-atomic separation. For weak pin- 
ning, we can still use the rigid approximation for skyrmions. 
In this case the interaction energy is just the self-energy of the 
skyrmion given by using J^xir) in Eq. ( [TT|) . We first obtain 
the structure of the skyrmion from Eq. ( |10| i and then calcu- 
late its self-energy with Eq. ([3]). We also do a full numerical 
relaxation by holding the spin at the center of the skyrmion 
unchanged in order to pin the skyrmion at a desired position. 
Both methods yield consistent results, as shown in Fig. |4] 
Several observations are as follows: 1) The length scale of 
the exponentially decaying force at a large distance is given 
by the size of defects Because the main contribution to 
the self-energy comes from the core of the skyrmion, the ex- 
ponential tail does not contribute significantly and the inter- 
action range is determined by ^j. 2) The force is maximized 
when the separation becomes close to the skyrmion radius, 
rd ~ Rc, and it drops when the skyrmion gets even closer to 
the defect and finally vanishes when the core coincides with 
the center of the defect. 3) Because the amplitude of the force 
is proportional to the strength of the defects, the force can 
be expressed as ~ Jdexp(-rd/^d) for large separations. 
4) The force is attractive for Jd < and it is repulsive for 
Jd > 0, so skyrmions prefer to stay in the Jd < region to 
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FIG. 5. (color online) (a) and (b) Real-space configuration of the 
skyrmion lattice in the absence of temperature gradient (a) and in the 
presence of temperature gradient (b). For clarity, only part of the 
configuration is shown, (c) Bragg peaks (green dots and lines) of the 
skyrmion configurations with a temperature gradient. The triangu- 
lar skyrmion lattice is rotated and distorted. The Bragg peaks for a 
perfect triangular lattice without a temperature gradient (orange lines 
and dots) are shown for comparison. 



5 



dient. flT\ The rotation was explained in Refs.l23]andl24lbv 
using a continuum model. Here we show that the rotation can 
also be explained by using our particle model in Eq. Q. 

Fjs = Frf = in the crystal phase without defects. The Hall 
angle of the skyrmion trajectory is tanOn = Vy/v^ = -a/rj 
when the current is along the x direction. The damping coef- 
ficient, a, has two contributions: the Gilbert damping and the 
dissipation due to the electric field induced by the skyrmion 
motion ll22l . The latter contribution is dominant. Thus, 
a ~ (T, where cr(r) is the temperature-dependent conductiv- 
ity. The temperature gradient leads to a spatial variation of the 
Hall angle that exerts a finite torque on the skyrmion lattice. In 
the absence of pinning, the lattice keeps rotating with a finite 
angular velocity. However, pinning is always present real sys- 
tems due to the underlying atomic crystal structure that favors 
one particular orientation of the skyrmion crystal. An addi- 
tional pinning arises from the geometric confinement of finite 
samples. The competition between torque and crystal pinning 
thus yields to a stationary state in which the skyrmion lattice 
is rotated by a finite angle. The lattice keeps rotating with fi- 
nite angular velocity for a sufficiently large torque induced by 
large enough currents and/or temperature gradients. 

We perform numerical simulations of Eq. (|5]) by modeling 
the temperature gradient with a(x) - 0.6 - 0.5(2x - Lx)^ IL^, 
where Lj is the length along the x direction. The current is 
also parallel to the x direction and we use a simulation box 
with aspect ratio : Ly - 2 : V3 and periodic boundary 
conditions. These boundary conditions favor a particular ori- 
entation of the principal axis of the skyrmion lattice, which 
is parallel to the x direction. In the stationary state, we find 
that the skyrmion lattice is rotated by a finite angle relative to 
the case of zero temperature gradient ao = 0.1 [see Fig. |5). 
Thus the experimental observation can be explained with the 
particle model in Eq. (|5]). 



B. Creep motion of skyrmions 

The skyrmions can easily leave the pinning potential either 
by quantum or/and thermal fluctuations because pinning is 
weak. This phenomenon leads to creep motion. We consider 
the dynamics of a single skyrmion in a pinning potential U(r). 
Because a <K 1 for real materials, such as MnSi, we will ne- 
glect the damping for simplicity. The quantum creep rate for 
Eq. (|5]l was calculated in Refs. [T5]and|25]for superconducting 
vortices. To be specific, we will consider a pinning potential 
per unit length 



/ x^ x^ 



The action of Eq. (|6]) becomes 



dT 



An 
7 



-ydjX ■ 



^2 



X 
i2 



X 



(12) 



(13) 



in the imaginary time representation t — > -It. Here y 
plays the role of momentum; thus, if the potential is sep- 
arable, i.e. U{x,y) - U\{x) + U2iy), the y dependent po- 
tential is not inverted in the imaginary time representation. 



Equation ( pj) ) is the same as the one for a particle with 
mass m - A^(2Ud) moving in a one dimensional potential 
U(x)/Ud - x^/A^ - x^/A^. The skyrmion does not have an in- 
trinsic mass according to Eq. (|5]). However, it gains an extrin- 
sic mass in the presence of a pinning potential. The quantum 
rate, Tg ~ exp{-Sg/fi), with 



n 



32nA^d 
15yn ' 



(14) 



is independent of the height of the pinning potential, but it 
depends on the width, which is a generic feature of quantum 
tunneling. 

We now consider the escape rate due to thermal fluctua- 
tions. For this purpose we add a noise force F„ in Eq. ^ 
which satisfies 

Ana 

(F„) = 0, {F„(t)F„{t'))^2kBT 6(t-t'), (15) 

m 

according to the fluctuation-dissipation theorem. The thermal 
rate, Yj - Q.exp(-AU/kBT), is dominated by the exponential 
factor, exp(-AU/kBT), where AU - 4Udd/21 is the height 
of the pinning potential. This factor reflects the Boltzmann 
distribution of the skyrmion in the potential U, and is thus 
independent of the dynamics (see Appendix C). In contrast, 
the attempt frequency, Q., does depend on the dynamics and 
the Magnus force. 

Thermal escape becomes dominant at high temperatures, 
while quantum creep is dominant in the low temperature re- 
gion. The crossover temperature between quantum and ther- 
mal tunneling is keT* = AUh/Sg = 5ynUd/(72nA^), which 
depends on the ratio of the width of the pinning potential and 
its height. T* can be increased for properly engineered pin- 
ning potentials. Since the skyrmion carries a magnetization 
that is opposite to the ferromagnetic background, the creep 
motion manifests itself in experiments as a decay of the oppo- 
site magnetization. Thus, like in the case of superconducting 
vortices, the rate can be obtained by measuring the time de- 
pendence of the magnetization. 



C. Comparison between continuum and particle models 

To validate the particle model, we perform numerical sim- 
ulations with both the particle and continuum models. We 
calculate the velocity of skyrmions as a function of the driv- 
ing force when defects are present. The defects are modeled 
as in Eq. ( [TT] i. Skyrmions are pinned in a low driving cur- 
rent, and they depin from the defects when the Lorentz force 
is high enough. The particle model yields results that are in 
reasonable agreement with the continuum model. Near depin- 
ning, the electric field behaves as E^j ~ (J - Jcf'', where Jc 
is the depinning current. From the numerical data we obtain 
Px ~ 1.55 and ySj ~ 1.93. The exponent /3f^ > I indicates that 
the depinning is plastic, i.e., some skyrmions escape from the 
pinning centers, while the others remain pinned. Eventually, 
all skyrmions become depinned when the current is further 
increased. 
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FIG. 6. (color online) Comparison between the continuum and parti- 
cle models. The lines are obtained with the continuum model while 
the symbols are obtained with the particle model. We convert the 
velocity of skyrmions in the particle model to an electric field by 
multiplying by an appropriate scaling factor The inset is the scaling 
of the electric field near the depinning transition with = 0.007. 
The lines in the inset are power-law fits. 



It is also interesting to discuss the effect of the Magnus 
force on the pinning of skyrmions. The Magnus force domi- 
nates over the dissipative force, Fm » 47rav/(riy), for a/rj <g: 
1. When a skyrmion moves around a pinning center or an ob- 
stacle, it is easily scattered with a velocity perpendicular to the 
pinning force or repulsive force. Thus, the skyrmion avoids 
passing through the pinning center and its influence is mini- 
mized, as shown in Fig.|7](a) and (b). When the dissipative 
force is dominant, a/rj » 1, the skyrmion has to pass through 
the pinning center so the pinning becomes very strong, as it 
shown in Fig. |7jc). We performed numerical simulations for 
several a/i] ratios [see Fig. |7](d)]. The depinning current is 
weaker for smaller values of a/j], i.e., when the Magnus force 
is dominant. 



IV. DISCUSSION 

The equation of motion for the center of mass of a rigid 
skyrmion lattice has been derived in Ref. 24 by using an ap- 
proach proposed by Thiele [26J. This equation is similar to 
Eq. (|5]l for F,,, = and Fd - 0. The equation of motion for 
the collective excitations in the skyrmion lattice was derived 
in Ref. 22, and it also shares a similar structure with Eq. (|5]l. 
The equation of motion presented in this work is more general 
because it also describes the dynamics of a single skyrmion 
and its interaction with other skyrmions and defects. 

The resulting equation of motion is similar to that of 
vortices in type II superconductors. However, the pinning 
of superconducting vortices is much stronger than that of 
skyrmions because of the following reasons. The Magnus 
force is negligibly small for vortices, except in the super clean 



region, which has not yet been realized experimentally. ifTSll 
In contrast, the Magnus force of skyrmions is stronger than the 
dissipative forces. Vortices have to pass through the pinning 
centers, which leads to large critical currents, while skyrmions 
can be easily deflected by the pinning center because of the 
dominant Magnus force. In addition, defects suppress super- 
conductivity over a range that is equal or larger than the su- 
perconducting coherence length (linear size of the normal vor- 
tex core). This matching of length scales makes the pinning 
rather strong. In contrast, the characteristic length of defects 
(inter-atomic spacing) is much smaller than the core size of 
the skyrmions. 

In conclusion, we have derived an effective particle model 
for skyrmions, which includes repulsive skyrmion-skyrmion 
interactions, interaction with defects and the role of the Mag- 
nus force. The model successfully describes the rotation of 
the skyrmion lattice in the presence of a temperature gradient 
and explains the small depinning thresholds that have been ex- 
perimentally observed. It also provides clear predictions for 
quantum and thermal creep. Finally the model has been vali- 
dated by direct comparisons resulting depinning and transport 
curves against the original continuum model. Our particle 
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FIG. 7. (color online) (a) and (b): Schematic view of a skyrmion 
passing through a pinning center (a) and obstacle (b) when the Mag- 
nus force is dominant. When the Magnus force is dominant over the 
dissipative force, the skyrmion is deflected by the pinning centers or 
obstacles, (c) Same as (a) and (b) except that the dissipative force is 
dominant. The skyrmion has to overcome the pinning site or obsta- 
cle by passing through it. (d) Numerical results of current-velocity 
curves with different a/rj. The depinning current increases with a/rj. 
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FIG. 8. (color online) The frequency of the different indexed normal 
modes. Inset: the frequency of the lowest five modes at different /?. 



model offers a convenient and transparent theoretical frame- 
work for the future computational and analytical studies of 
skyrmions. 
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Note added - After completion of this work, we becomes 
aware of the work by Iwasaki et. a/.. ll27l where they also 
found using a continuum model that the pinning is weakened 
due to the presence of Magnus force. The dynamics of a 
skyrmion in the presence of a defect was considered in Ref. 



Appendix A: Internal mode of a skyrmion 

Here we study the normal modes associated with small fluc- 
tuations of the spin texture around a stationary skyrmion. In 
the polar coordinate (r, 0), the spins are described by 

n(r, 0) = (sin0cos(i^ - 0)r, sin0sin(^ - 0)0, cos0z), (Al) 

where f , and z are the unit vectors along the correspond- 
ing axises. We consider the fluctuations that only depend on 
r, which are expected to have lower energy. This approxima- 
tion excludes the mode corresponding to the translational mo- 
tion of the skyrmion. We consider the localized modes in the 
core region of the skyrmion, where spins are strongly canted 



from the direction of the applied magnetic field. We write 
Q{r) - do(r) + 6(r) and (p(r) - (p + njl + lp(r), where 0o(r) is 
the solution for a stationary unperturbed skyrmion. The weak 
perturbations 6(r) <K 1 and if(r) <K 1 are governed by 

d,0 + a sin 0od,(p = IjD sin OodrOoip 

J J 
+y-[(sm0Q + 2rcos 0odr6o)dr(p + rsm(0o)d^(fi], (A2) 
r 



ad,0 - sin 
2£)sin(20o) ycos(20o) 



r 



- H„ cos 01 



(A3) 



Equations ( |A2[ ) and ( |A3| l are not applicable in description of 
the extended modes (spin wave excitations) where the expan- 
sion around 00 is not valid because 0o <ii 1. For the continuum 
spin wave excitation, the structure of skyrmion becomes irrel- 
evant, and the dispersion is given hy co - /exk^ + yHg with- 
out damping a - 0. To obtain the spectrum of the localized 
modes, first we obtain 0{){r) by solving Eq. ( [TO] i numerically, 
and then we solve Eqs. ( |A2[ i and ( |A3[ ) with the boundary con- 
dition 0(r = 0) = ifi(r = 0) = 0(r - oo) - ip(r = oo) = 0. 
We neglect the damping and take a - 0. The results for the 
spectrum and the lowest eigen modes of the normal modes are 
shown in Figs. [8]and[9j where /3 = 2HaJex/D^- The normal 
modes are always gapped for a nonzero D. The gap increases 
with magnetic fields and is larger than the spin-wave gap for 
a large D. However, for a sufficiently small D, the gap of the 
normal modes becomes small and vanishes at D = while the 
spin-wave gap is independent of D. The lowest normal modes 
lie below the spin wave excitations in this case. The eigen 
modes can be describe by damped oscillatory functions. The 
period of oscillation decreases and the frequency of the modes 
increases. Using typical parameters for MnSi \22\, /3 ~ 1, the 
spin wave gap is 10 GHz and the normal modes have even 
higher frequency. Thus the structure of the skyrmion is pro- 
tected by the existence of the gap in the normal modes. 
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FIG. 9. (color online) The eigen mode for the lowest three modes at 
/?= 1.2. 
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Appendix B: Numerical details 

In simulations, we introduce dimensionless units in Eqs. (|3]l 
and Q. Length is in unit of Jex/D; energy is in unit of Jg^/D; 
magnetic field is in unit of jj^^; time is in unit of J^^lijD^); 
current is in unit of IDejti. We use the periodic boundary con- 
dition in both directions. To find the ground state, we anneal 
the system by adding a Gaussian noisy magnetic field along 
the z direction in Hgff. Equation Q is solved by an explicit 
numerical scheme developed in Ref. 29 The current is along 
the X direction. In calculations of the results in Fig. [3] we use 
a simulation box of size x Ly = 30 x 10. The system is dis- 
cretized with a grid size of 0.2. In calculations of the results 
in Fig. |6] the defects are modeled by Eq. ( [TT| i with Jd - 1.0 
and - 1.0, where the interaction between skyrmions and 
defects is repulsive. The Nj - 500 defects are randomly dis- 
tributed in a simulation box of size X Lj, = 100 x 100. The 
I-V curves are obtained by averaging over 20 realizations of 
random defects. 

In the particle-level simulation, we take the interaction 
between skyrmions as F,, = FsCiKiir^/^)?^ and the repul- 
sive interaction between skyrmions and defects as - 
FciQQxp{-rdl^d)rd- Here r^/ is a unit vector along rd- In di- 
mensionless units, the equation of motion becomes 



"17 



In simulation = 2^ = 2.0, Fso = 1.0 and Fda - 0.6. The 
number of skyrmions and the number of defects Nd aie N - 
Nd - 225. The simulation box is L, x L,, = 45 x 39. We use 
the second-order Runge-Kutta method to integrate Eq. ( |B1| | 
with a time step At - 0.05. 



Appendix C: Thermal activation of a skyrmion over a barrier 

Here we calculate the thermal activation rate for a skyrmion 
in a metastable potential U{x,y). The equation of motion for 
the skyrmion in dimensionless units is 

i+[zx\-VU]^a\, (CI) 

with a Gaussian noisy force (/^(r,f)) =Oand 

{Mr,t)f,,(r',t')) = 2aT6ir-r')6it-t')6^,,,. (C2) 



the initial position ro at to is 

/?(r',f'|ro,fo) = j D[r,f]exp(^j dtJ^, (C3) 



--drU -H -Z X V 

a a 



(C4) 



with T - 2T /a. It is more convenient to use the Hamiltonian 
description 



Aa^ [pI + p]) d,U [ap, + p,) dyU [ap, - p,) 



2(Qf2 + 1) 



+ 1 



0-2 + 1 



Using the nonequilibrium path integral approach 113011 . the 
probability of finding the skyrmion at r' at t' starting from 



where the conjugate momentum p is defined as p ^ = dX./dvx 
and py - dX./dvy. For a weak noise f <s; 1, the dominant con- 
tribution to the path integral are those trajectories governed 
by the standard Hamiltonian dynamics. For the Hamiltonian 
dynamics, "7/ is conserved. Initially for the skyrmion at the 
well r„, the system has 'H - since p^ = py = 0. We then 
look for the trajectories with 'H - Q and with minimal action 
St - - J dtX.. One obvious solution for 9i - Ois p^ - Py - 
with 5 r = 0. This is a non-fluctuating trajectory, which does 
not contribute to the thermal activation of skyrmions. There is 
another trajectory with a minimal S t 



P, = -i=r^,U, (C5) 
Ta 

with fu - X, y. The corresponding equation of motion is 

[z X V - Vt/] = -a\. (C6) 

Compared with Eq. ( |Cl| l, the sign of damping is changed. 
Thus Eq. ( |C6| l describes the motion of a skyrmion in a po- 
tential with a negative damping, which forces the skyrmion to 
leave the well, and contributes to the thermal activation. The 
action S j for this trajectory is 

r 2 r'" Au 

S^^ j dtpy^-^ j [VU-dr]^—, (C7) 

where the integration is from the well r„, to the barrier r;, of 
the potential U, and AU - Uivt) - C/(r,„) is the height of 
the barrier. The action does not depend on the Magnus force. 
The reasons are as follows: first the dynamics of skyrmions is 
irrelevant for the thermal activation of skyrmions over the bar- 
rier. The probability distribution of skyrmion in the potential 
only depends on the potential energy. Secondly, the Magnus 
force does not produce work when skyrmions move. 
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